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ABSTRACT 

Aims. This work presents a high-precision variability survey in the field of the old, super metal-rich open cluster NGC 6791. 
Methods. The data sample consists of more than 75,000 high-precision CCD time series measurements in the V band obtained mainly 
at the Canada-France-Hawaii Telescope, with additional data from S. Pedro Martir and Loiano observatories, over a time span of ten 
nights. The field covers an area of 42 x 28 arcmin 2 . 

Results. We have discovered 260 new variables and re-determined periods and amplitudes of 70 known variable stars. By means of 
a photometric evaluation of the membership in NGC 679 1 , and a preliminary membership based on the proper motions, we give a 
full description of the variable content of the cluster and surrounding field in the range 16 < V < 23.5. Accurate periods can be given 
for the variables with P <4.0 d, while for ones with longer periods the limited time-baseline hampered precise determinations. We 
categorized the entire sample as follows: 6 pulsating, 3 irregular, 3 cataclysmic, 89 rotational variables and 61 eclipsing systems; 
moreover, we detected 168 candidate variables for which we cannot give a variability class since their periods are much longer than 
our time baseline. 

Conclusions. On the basis of photometric considerations, and of the positions of the stars with respect to the center of the cluster, 
we inferred that 11 new variable stars are likely members of the cluster, for 22 stars the membership is doubtful and 137 are likely 
non-members. We also detected an outburst of about 3 mag in the light curve of a very faint blue star belonging to the cluster and we 
suggest that this star could be a new U Gem (dwarf nova) cataclysmic variable. 

Key words. Stars: starspots - Stars: statistics - Stars: variables: general - binaries: eclipsing - novae, cataclysmic variables - open 
clusters and associations: individual: NGC 6791 



1. Introduction 

The photometric precision achieved by several ongoing transit- 
ing planet searches allows us to extend the census of variable 
stars down to very low amplitudes and faint magnitudes in se- 
lected sky regions. Variable stars are an important source of as- 
trophysical information: from observations of them we are able 
to test several theories (e.g., evolutionary and pulsational mod- 
els). Since all stars in an open cluster have essentially the same 
age, chemical composition and distance, the study of variables 
which are cluster members can put more severe constraints on 
the physical parameters. Comparisons can also be made between 
the variable stars of the cluster and those of the surrounding field. 

In this paper, we present the study and the classification of 
260 new variable stars that we found in the field of the open 
cluster NGC 6791, while for 70 already known variables we 
compare our results with the previous ones. NGC 6791 (a- 19 h 
2Q m 53 s ; S= +37° 46' 18"), is a rich and well studied open clus- 
ter. It is thought to be the one of the oldest and probably the 
most metal-rich cluster known in our Galaxy. Its age is estimated 
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to be about 8-9 Gyrs (Carr aro et al., [20061 Kinget al., [2005 ; 
Chaboyer, G reen, & Liebert, [19991 Stetson et al., 2003) Kaluzny 
& Rucinski, |19951 l; however, the white dwarf cooling sequence 
indicates a different value, i.e., ~2.4 Gyr (Bedin et al., [2005). 
The most recent estimates of its metallicity are [Fe/H]=+0.39 
(Carrara et al.. 120061 [Fe/H]=+0.47 ( Gratton et al.. l2006l l. and 
[Fe/H]=+0.45 (Anthony-Twarog et al.. 120061 1. In this work we 
adopt for NGC 6791 a distance modulus (m-M) v = 13.3 5+0.20 
mag and a reddening E(B - V)=0.09 mag (Carraro et al., 2006). 
The cluster is thus located at about 4. 1 kpc from the Sun. 

Because of its extreme characteristics, NGC 6791 has been 
the target of many surveys (see Table[T]for a list of publications). 
Taking into account the fact that in four cases the same stars have 
two identification numbers (V15=B7, V56=V96, V76=V85 and 
V77sV88) and counting also the stars B4 and B8, the total num- 
ber of known variable stars in the field of NGC 6791 to date 
was 123 (plus 7 suspected variables, proposed by Hartman et al., 
120051 

In Sect. |2]below we describe our observations, in Sect. |3]we 
give details about the methods we employed in the search for 
variable stars, which are themselves presented in Sect. |3.2| In 
Sect. [4] we describe the properties of the variable stars, focusing 
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Table 1. Previous variable star searches in NGC 6791. 



Authors 


Nr. of variables 


IDs 


Notes 


Kaluzny & Rucinski (1 99311 (KR93) 


17 


V1-V17 


V15=B7 


Rucinski, Kaluzny & Hilditch i 1996b (RK96) 


5 


V18-V21 and B8 




Morheiska et al (200? 1 (MO?*! 


47 


V22-V67 and B4 


i_> t 1 w cis ijic-v luuuiy ciiciivjg ucu uy 

Kaluzny & Udalski (| 1992b. as 
a blue star, but not as variable. 


Mochejska et al. j2003t (M03) 


7 


V68-V74 




Kaluzny j2003l (K03) 


4 


V75-V78 




Bruntt et al. {US} (B03) 


19 


V79-V100 


V85=V76; V56=V96; V77=V88 


Mochejska et al. (2005 1 (M05) 


14 


V101-V114 




Hartman et al. ( 2005 ) 


10 


V115-V124 


Plus 7 suspected variables 



our attention on probable cluster members and some additional 
peculiar cases. The entire catalogue of variable stars is reported 
in an Appendix. 

2. Observations and data reduction 

We surveyed NGC 6791 to detect the transits of extrasolar plan- 
ets (Montalto et al., 12007b . The campaign covered 10 consec- 
utive nights (from July 4, 2002 to July 13, 2002) and it was 
characterized by the continuous monitoring of the target on each 
clear night. Therefore, in addition to the planetary transit search, 
we could get access to the full variability content at P <4.0 d, 
both for the cluster and the surrounding field. Three telescopes 
were used: 

1. The Canada-France-Hawaii Telescope (CFHT) in Hawaii 
equipped with the CFHT 12k detector, composed of 12 CCDs 
of 4128 x 2048 pixels and covering a field of about 0.32 
deg 2 . Owing to the large number of bad columns, data from 
chip 6 could not be used, so we could get data over a 
0.29 deg 2 field; 

2. The San San Pedro Martir (SPM) 2.1-m telescope equipped 
with the Thomson 2k detector and covering a field of about 
6x6 arcmin 2 ; 

3. The Loiano 1.5-m telescope equipped with BFOSC + 
the EEV 1300xl348B detector and covering a field of 
11.5 x 11.5 arcmin 2 . 

Table|2]gives details about the length of the observing nights 
while Figure[T]shows the field of the CFHT survey and the edges 
of the Loiano and SPM surveys. The coordinates of the edges of 
our fields are also listed in Tab. [2] The field of the SPM obser- 
vations is entirely included within chip 9 and the field of the 
Loiano observations partially covers chips 2, 3, 4, 8, 9 and 10 of 
CFHT (see Fig. [TJ. The luminosities of the new variables range 
from V=23.2 mag to V ~ 16 mag (near 1 mag above the turn- 
off); brighter stars are saturated. The calibration of the CFHT, 
Loiano and SPM data have been performed by using the Kaluzny 
& Rucinski photometry ( 1995). More details on the data reduc- 
tion procedure can be found in Montalto et al. ([2007). 

3. The identification of variable stars 

The intensive monitoring of NGC 6791 allowed us to obtain 
tens of thousands of photometric time series for stars located in, 
close to, and far away from the cluster center. We have analysed 
73331, 6055, and 2152 light curves obtained from the CFHT, 
Loiano and SPM telescopes respectively. The CFHT, Loiano and 
SPM time series are composed of about 250, 60 and 170 data- 
points, respectively. The observations, intended to detect photo- 
metric transits, were performed in the V band only. 



Table 2. The observation log for each night and the limits of the 
field of view at the 3 different observatories. 





Loiano 


SPM 


CFHT 


Night 


*-start tend 


t-start 


tend 


^-start 


tend 




(HJD-2452400) 


(HJD-2452400) 


(HJD-2452400) 


V 








59.82 


59.96 


2 nd 








60.83 


61.02 


yd 


62.48 62.63 


61.68 


61.95 


61.94 


62.07 


4* 


63.41 63.63 


62.68 


62.96 


62.76 


63.07 


5"' 


64.38 64.64 


63.68 


63.96 


63.88 


64.10 


6"' 


65.39 65.62 


64.69 


64.98 


64.77 


65.11 


nth 




65.70 


65.82 






w 




66.69 


66.97 






giii 




67.67 


67.98 


67.77 


68.10 


10"' 




68.69 


68.87 


68.76 


69.11 


&min 


19" 20™ 25 s .8 


19" 20" 


3ft3 


19" 19" 


23 s .7 


™max 


19" 21'" 30:4 


19' 21" 


lff.4 


19" 22" 


58*0 


&min 


37° 41'2276 


37° 43'15'.'4 


37° 36'6'.'7 


^max 


37° 53'42'.'7 


37° 50'3'.'l 


38° 4'21'.'2 



3.1. The search for variable candidates 

To search for variable stars, we calculated the best "sinusoid plus 
constant" fit for all light curves (Vanicek, I1971I Ferraz-Mello, 
1 1 98 lb - We evaluated the goodness of the fit by calculating pa- 
rameters related to the reduction of the initial variance obtained 
by introducing the periodic term. These parameters are the re- 
duction factor (Vanicek, 119711 1 and the coefficient of spectral 
correlation S (v) (Ferraz-Mello, 1 198 11 1. 

Owing to the huge number of light curves, we need one or 
more parameters to discover the variability. Toward this goal, we 
considered the parameter r defined as r — log I0 S m(lx , where S max 
is the maximum value of S (i.e., the one corresponding to the fre- 
quency of the best-fit sinusoid in the Ferraz-Mello method). If a 
star does not show variability the introduction of a sinusoid does 
not improve the fit and then S max is close to zero (no variance 
reduction) and r « 0; on the other hand, a sine-shaped variabil- 
ity strongly reduces the variance (S close to 1) and hence r = 0. 
The purpose was to use the r parameter as a tracer of variability 
for short-period (i.e., intranight) variability. 

To search for long-period variability, we introduced a sec- 
ond parameter, more sensitive to the night-to-night variations. 
We calculated the mean magnitude V, and the standard deviation 
cr, on each night, and after that we calculated the parameter s 
defined as: 



where AV is the peak-to-peak difference and cr is the mean of 
the cr,- over all nights. 
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Fig. 1. Field of view (42 x 28 arcmin 2 ) of the CFHT image. The chips are numbered in increasing order from left to right from chip 1 (top left) 
to chip 12 (bottom right); stars of chip 6 are not plotted since we found it impossible to derive accurate photometry. Dashed and solid lines are the 
edges of the SPM and Loiano fields, respectively. Variable stars are also plotted: blue dots are pulsating variables, purple dots are irregular and 
cataclysmic variables, green dots and cyan dots are eclipsing systems (EA/EB-Type and EW-Type, respectively). Finally, red and yellow dots are 
rotational and long-period variables, respectively. 



To test the capability of the r and s parameters to detect vari- 
able stars, we prepared a sample containing two types of light 
curves: 7722 artificial constant light curves (see Montalto et al., 
120071 for details) and 70 light curves of already known variable 
stars which are included in our CFHT field. In Fig. [2] we plot 
r vs. s for the light curves of constant stars (small points) and 
of variable stars (large points). The variable stars are substan- 
tially apart from the constant stars and most have r > — 1. The 
variable stars with r < — 1 and superposed on constant stars are 
mostly EA-type stars or irregular stars (e.g., cataclysmic vari- 
ables). Among variables (i.e., large dots in Fig.[2]i, the stars with 
small s have short periods (P < 0.50 d), while stars with large 
s have long periods. Therefore, we can conclude that the combi- 
nation of the r and s parameters is a good tracer of variability. 

To detect the variable stars in our sample of ~82,000 light 
curves we first selected in an automatic way all the stars with 
r > -2.0, according to the test described above. We thereby re- 
duced the huge initial sample to ~6,500 stars. After the calcula- 
tion of the amplitude spectrum of their time series, we adopted as 
a second selection criterion a signal-to-noise ratio (S/N) greater 
than 4.0 around the highest peak in the amplitude spectrum. This 
procedure allowed us to reduce our sample to ~900 stars, i.e., 
1.1% of the whole initial sample. Further checks have been made 
by examining the light curves of a random sample of stars with 
r < -2.0, large s and 3.5<S/N<4.0, but we did not find any ad- 
ditional variables. 

Our approach allowed us to detected hundreds of stars show- 
ing peaks in their power spectra at / =1.00 cd~', at / <0.05 cd~\ 
or at /=0.6 cd~'. The first two spurious periodicities are com- 



mon and can be ascribed to small misalignments in the mean 
magnitudes from one night to the next or recurrent drifts (caused 
by small color effects, for example) in the intranight light curves. 
We suggest that the latter one is probably a photometric artefact 
occurring in some particular cases of blended stars, or stars close 
to CCD edges, or bad pixels. They have been considered as not 
reliable enough to infer a physical light variability. In our opin- 
ion, only the combination of automatic procedures and visual 
inspection allowed us to identify the three classes OP=1.00 d, 
P=\.6 d, P »10 d) of spurious variables in the huge number of 
~82,000 light curves. In particular, we note that the identification 
of the whole sample of eclipsing binaries has been confirmed by 
the application of the box fitting technique (BLS, Kovacs et al. 
2002 ), used by Montalto et al. d20071 > to detect planetary transits. 

At the end of the variable star identification, we were left 
with 330 cases to be characterized. Since we rejected about 2/3 
of the sample selected by means of the r, s parameters, we are 
confident we have not applied overly strict constraints in the can- 
didate selection. 

3.2. The cross check with previous surveys of NGC 6791 

When comparing our field of view with those of other surveys, 
we found that 8 1 known variable stars are included. The CFHT 
survey failed to detect 45 known variable stars: seventeen stars 
(V22, V24, V26, V28, V30, V35, V36, V47, V50, V57, V61, 
V63, V64, V102, V103, V104, V105) are outside the CFHT 
field of view; 4 stars (V71, VI 06, VI 13 and VI 20) lie between 
two chips; 23 stars (VI, V6, V13, V19, V33, V39, V45, V49, 
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Fig. 2. Grey points: r parameter vs. s parameter for constant light 
curves. Black dots: parameters derived from our light curves for the 
variable stars previously detected in our field. 



V54, V56=V96, V65, V66, V67, V69, V70, V72, V73, V74, 
V77=V88, V78, V81, V97 and VI 12) are saturated; and the V76 
data are useless. 

Among the 81 known variable stars that we have observed, 
not all of them display variability in our sample: 4 stars (V10, 
VI 8, V21, V32) are previously classified as long -period de- 
tached eclipsing variables and we did not observed eclipses. We 
are not able to confirm the period of 15.24 days for V68 (M03), 
likely because of our shorter time baseline and the small ampli- 
tude of this variable (about 0.003 mag in V-band, M03). Finally, 
we cannot confirm the variability of six stars (V20, V79, V84, 
V98, V99, VI 16) and of the seven suspected variables found by 
H05, since our data do not show any trace of variability. 

Among the sample of the stars missing from the CFHT field, 
we identified 22 stars in the Loiano and SPM data sets (V6, 
V13, V19, V20, V33, V45, V54, V56=V96, V65, V66, V67, 
V70, V71, V73, V74, V76=V85, V77, V78, V81, V97, V106 
and VI 13). However, owing to the smaller signal-to-noise ra- 
tio (S/N), the small number of datapoints and (in the case of the 
Loiano data) the limited survey time, we could only confirm the 
variability of stars V56=V96, V66 and V76=V85. 

Throughout this paper we use the existing names for the 
already known variables; to identify the new ones discovered 
in our survey we used the five-digit number assigned by the 
DAOPHOT package followed by the number of the chip which 
the star belongs to. Accurate astrometry is provided to identify 
the stars on the sky. Moreover, all light curves of the variables 
will be available on CDS. 



4. The variable star content of NGC 6791 and its 
surrounding field 

The CFHT measurements are quite precise, thus the light curves 
are generally very well defined for P <4 d. On the other hand, the 
periods and the shapes are uncertain for P >4 d, since the obser- 



Fig. 3. Proper-motion vector-point diagram for the inner region of 
NGC 6791 (from Bedin et al., 20061. The circle (centered on the abso- 
lute proper motion of the cluster) represents a safe limit corresponding 
to 0.5 mas/yr. Triangles represent non-members, points with error bars 
are the new variables 06289_9, 04803_9 and 0983 1_9 



vations only covered 2.5 cycles or less. We refer to Montalto et 
al. (J2007) for a full description of the photometric errors. In or- 
der to evaluate the precision in the study of the variable stars, we 
calculated the standard deviations of the Fourier least-squares 
fits (truncated at the last significant term for the given star) for 
the 138 light curves having very good phase coverage. The pre- 
cision was found to be better than 0.010 mag in 73 cases (53%), 
and better than 0.020 mag in a total of 122 cases (88%), as ex- 
pected for stars ranging from V ~16.0 to V ~22.5. The discus- 
sion based is mostly on the CFHT data, which are by far the 
most numerous, precise and homogeneous; however, for some 
variables we have used data from Loiano and SPM in a very 
profitable way. As an example, only the longitude spread of the 
three observatories allowed us to derive the periods of the eclips- 
ing binaries 00645_10, V107, V12, V109 and of the rotational 
variable 03079_9. 

To proceed in the definition of the variable star content of 
NGC 6791 and its surrounding field, we calculated the power 
spectra of the data for all the 330 candidate variables by using the 
least-squares iterative sine -wave search (Vanicek. 119711 1 and the 
Phase Dispersion Minimization (Stellingwerf, 1978) methods. 
Differences have been examined and resolved. The separation 
into different classes of variable stars has been made on the basis 
of the light curve parameters (period, amplitude, Fourier coeffi- 
cients) and standard photometric values (V, B - V, V - I), when 
available. The period estimates have been refined by means of 
a least-squares procedure (MTRAP, Carpino ll987t and appro- 
priate error bars have also been calculated. At the end of the 
process we get six pulsating stars with P < 0.6 d, three irregu- 
lar variables, three cataclysmic variables (CVs), 31 detached or 
semi-detached eclipsing binaries, 29 contact binaries, 90 rota- 
tional variables, 167 stars showing clear night-to-night variabil- 
ity on timescales too long for periods to be determined over our 
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9.2-d baseline. We adopt preliminary membership probabilities 
based on proper motion measurements kindly provided to us by 
K. Cudworth (private communication) for 35 stars. Moreover, 
for three new variable stars we adopted membership probabili- 
ties based on proper motions performed by Bedin et al. (2006) 
(hereafter B06, see Figure [3]). 

For the other stars, we consider their position in color- 
magnitude diagrams (CMDs), and their distance from the cen- 
ter of the cluster to infer whether they belong to the cluster (for 
EW-Type stars we also utilize the P-L-C relation of Rucinski 
(2003 ). Toward this end, we plotted the radial distribution of all 
stars in Fig. [4] We see that at a distance of ~10' from the center 
of the cluster, the stellar density becomes near constant (about 
21 star/arcmin 2 ). Thus, we adopt the value of 10' as the external 
limit of the cluster and we consider "likely non-members" the 
variables located farther from the cluster center. 



100 



1 



a 

Z 



o. 






5 10 15 

distance from the center (arcmin) 

Fig. 4. Stellar density p (number of stars per square arcminute) as a 
function of the distance from the center. The straight line represents the 
mean stellar density at distances greater than 10'. 



Moreover, there are several variables whose light curves 
are very similar to those of Cepheid variables; the Fourier de- 
composition of some light curves (in particular the large am- 
plitude ones, i.e., 00913_5, 01659_8, V46 and 01431J0, but 
also 01606_11, 02285_10, 00122.4 and 03056_3) yields param- 
eters typical for Cepheid light curves. However, most of these 
variables are quite faint and the Period-Luminosity relation for 
Cepheids (Tammann et al., (2003 ) yields distances in the range 
39-171 Kpc. It is difficult to say whether these stars are nearby 
rotational variables (see below) in the Milky Way or very dis- 
tant pulsating variables. For our present purposes, these stars 
have been included among the rotational variables listed in the 
Appendix. 

The puzzling nature of all these apparently distant stars (i.e., 
the Cepheid-like ones, the RR Lyr and the faint SX Phe vari- 
ables discussed above) deserves further investigation by means 
of spectroscopic and/or kinematic data. 

4.2. Irregular variables 

Table [3] also lists three irregular variables: these stars lie on the 
middle Main Sequence and are all located less than 3' from the 
cluster center; thus we suggest that they belong to the cluster. 
V92 and V83 were previously defined as "periodic variables" 
by B03. Indeed, we noticed fast variability in our light curves 
(Fig. [6}, but, more noticeably, the mean magnitude is also chang- 
ing from night to night. The long periods given by M03 are not 
able to explain either the short- or the longer- timesc ale variabil- 
ity; actually, we could not detect any periodic term. We also de- 
tected no trace of periodicity in V93 (Fig. [6}; we suspect that 
the periods given by M05 and B03 are spurious, since they are 
close to 1.0 d (0.99 and 0.94 d, respectively) and they could be 
produced by the irregular fluctuations. 

We can conjecture that these variables are eruptive vari- 
ables observed in a quiescent phase, in which rapid and/or slow 
changes with smaller amplitude can be observed; they resemble 
the case of V15 (see Sect. [43}. We have no reliable indications 
about the membership probabilities. 



4.1. Pulsating variables 

The main characteristics of our variables are listed in Tab. [3] and 
their light curves are shown in Fig. [5] The classification as High- 
Amplitude Delta Set (HADS), SX Phe, RRc or RRab stars is 
based on the parameters of the Fourier decomposition (Poretti, 
2001). In all cases, the 021 Fourier parameters are on the progres- 
sions described by the different classes. We note that our period 
for V123 is quite different from that given by H05 (0.107 d). 
Error bars on the periods are in the range 1-6-10 -5 d. 

Both RR Lyr variables are too faint to belong to NGC 6791. 
Since they have V=17.21 (03653_3) and V=18.28 (00345_1), 
their distance moduli greatly exceed that of the cluster. 

This is also true for the very faint and short-period stars 
0031 IJ (7=23.17) and 00224J0 (V=21.72); therefore, it is 
more likely that they are Pop. II stars. On the other hand, using 
the P - L relation given by McNamara (2000), we get distance 
moduli of 14.50 and 13.78, respectively, for V123 and 01497J2. 
These distance moduli and the distance from the cluster center 
(12' and 22', respectively) suggest that they do not belong to the 
cluster, though they are not very far from it. Therefore, they are 
probably Pop. I stars and hence High Amplitude 5 Scuti stars. 



4.3. Cataclysmic variables 

As regards V15: M03 and M05 detected variability over the 
range of 3 mag and observed outbursts of about 0.5-1.0 mag; 
from our side, we could see a 0. 15-mag variability in our light 
curves (Fig. Eft, corresponding to the quiescent phase. V15 is 
very probably a NGC 6971 member, since the Cudworth proper- 
motion membership probability is very high (98%). 

Both the position of the faint blue star 06289-9 in the two- 
colour diagram and the shape of its light curve (Figure |7]i 
strongly suggest that this star could be a new cataclysmic vari- 
able (U Gem-type, dwarf nova). Moreover, we know that this 
object is a cluster member (see Figure [3}. The star shows an 
outburst of about 3 mag and, though we did not observe the en- 
tire brightening, we would highlight that the magnitude was still 
increasing on the first night; thus we are able to say that the max- 
imum brightness was reached immediately after. 

We can estimate the orbital period, P or b, and the recur- 
rence time, T„, from the decay time, r^-At/Am [days mag -1 ] 
and the amplitude, Am (Warner, 1995, equations 3.5, 3.1, re- 
spectively). Assuming for At and Am the values 3.33±0.50 d 
and 2.87±0.31 mag respectively, we find P ,./,=2.54±1.41 h and 
7'„ = 13.9+10.6 d. However, our light curve (Pig. |7j seems to rule 
out T„ values shorter than 8 d. 



6 



De Marchi et al.: Variable stars in the open cluster NGC 6791 and its surrounding field 



The variable B8 shows a large-amplitude light curve (Fig. |7]» 
over a quite short 7 d time span. The cataclysmic nature of B8 
has been confirmed spectroscopically by Kaluzny et al. (1997) 
who also notes that B8 exhibits red V - I colour while in a low 
state. 

Following the same procedure used for 06289_9 and assum- 
ing r d =l. 3+0.3 d mag' 1 for B8, we find P orfe =2.97±1.63 h and 
T n =\ 1.4+8.5 d. The T„ value is compatible with the 7 d period- 
icity (Fig. |7j. A membership probability is not available for B8. 
However, using the equations 3.3 and 3.4 after Warner dl9951 l. 
we obtain M Vmin -S. 06+0.68 mag and My mflJC =4.97±0.42 mag. 
In turn, these values give two estimates for the distance mod- 
ulus of B8, i.e., 13.82 + 0.68 mag and 14.20 ± 0.42 mag. We 
note that the first is in agreement with the distance modulus of 
the cluster. Kaluzny et al. d 19971) assumed that B8 belongs to the 
cluster, finding My„ M . v =5.2 mag and My m ,„=7.6, i.e., values very 
similar to ours. B8 is located at 4' from the center, and we can 
only conclude that the membership of this star is very probable. 
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Fig. 5. Light curves of pulsating variables. V123 and 01497_12 are 
probably HADS stars, 0031 1_7 and 00224.10 are SX Phe stars, 03653_3 
and 00345_1 are RR Lyr stars. 



4.4. Contact binaries 

The simplest cases of eclipsing systems are the contact bina- 
ries (also named W UMa systems); they show short periods and 
continuous variability and therefore can be easily recognized 
and classified. We detected 29 of these variable stars; they have 
P <0.40 days and very well defined light curves. The complete 
list and the light curves are reported in the Appendix. Tab. [4] 
lists the stars likely belonging to NGC 6791 (see above); their 
light curves are shown in Fig. [8] The very short periods and the 
secondary minima occurring at (f> as 0.5 indicate binaries with 
circular orbits, as is also the case for stars with small amplitudes 
(in 7 cases we have amplitudes less than 0.20 mag: V3, V4, V5, 
V8, V23, V40 and 01441.8). The average error bar on the period 
estimates is of the order of 4-5- 10~ 5 d. 

However, we note that the stellar surfaces are not homoge- 
neous since the maxima are often at different heights. Therefore, 
binarity and activity are probably combined here. In particular, 
the shape of the light curves of V4 (comparing RK96, M02 and 
our data) and V7 (comparing K93 and our data) have changed 
a lot; we suppose that stellar spots strongly modify the light 
curves. Proximity effects are also responsible for the large am- 
plitudes observed for 01434_3 (Fig.[|) last panel) and 00766_5. 
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We also found different periods for VI 18 (0.306321 d) and VI 24 
(0.320143 d) compared to H05. 

As for membership, the probabilities provided by Cudworth 
are 78%, 98% and 98% for V3, V4 and V5, respectively. Indeed, 
they are very close to the cluster center (4'5, 2! 1 and 2!8, respec- 
tively). 

We suggest that 01441_8, VI 18, V8, VI 17 and V7 are also 
contact binaries belonging to NGC 6791. To further verify this 
hint, we calculated their distance by using the P-L-C relation 
given by Rucinski (2003 ); they turn out to have distance moduli 
(13.28, 13.28, 13.48, 13.28 and 13.18, respectively) very similar 
to that of the cluster (13.35). 

Moreover, these stars are located at similar angular distances 
from the cluster center (6!2, 7'2, TA, T.2 and 6!3, respectively). 
Figure [9] shows how the distance modulus of the cluster is in 
better agreement with those of the stars we proposed as clus- 
ter members than with those of the previously known members. 
Their positions in the CMDs (Fig. 10 filled circles) are similar 



to those of stars in the sample with V <7.5 whose parallaxes 
have been determinated by HIPPARCOS (Rucinski, 120031 . We 
also note that most of the cluster members are near the turnoff 
point. 

4.5. Eclipsing variables 

In the cases of detached or semi-detached eclipsing binaries the 
classification and membership tasks are different from the case 
of contact binaries. Tab. [5] lists the systems for which we could 
determine periods; their light curves are shown in Fig. 1 1 We 
still have short-period cases where we can reconstruct the com- 
plete light curve, as for the classical examples of j3 Lyr variables 
(V29, 01558_5 and 00331_3). V9 is a more complicated /3 Lyr 
system in which spots produce maxima with different heights. 
Indeed, it has been classified as an RS CVn variable by M05 and 
B03; they also observed a "shift of the modulation wave" from 
1995 to 2002. 
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Table 3. Pulsating, irregular and cataclysmic variables. V is the minimum brightness for CVs and irregular, the mean brightness 
for pulsating variables. 7q is the time of maximum brightness for pulsating stars. Hereafter, the labels "k" and "s" indicate that the 
B - V color index is taken from Kaluzny & Rucinski d!995l or Stetson et al. (2003), respectively. 



Star 


Type 


Q '2000 


<52000 


V < B-V > 


<V - I > 


Ref. 


To 


Period 


Ampl. 










[mag] [mag] 


[mag] 




[HJD-2452400] 


[d] 


[mag] 










Pulsating variables 












V123 


HADS 


19.362064 


37.666034 


17.08 0.45 




k 


59.559 


0.06026 


0.14 


n 1/107 M 
U 14 7 / _1Z 






•27 Q 1 OA 1 Q 
J 1 .0 IZ'+IV 


10. UD 






jy.DZo 


u.u / zz / 


n ac\ 
u.^u 


00311-7 


SXPhe 


19.324628 


37.716768 


23.17 






59.605 


0.10443 


0.10 


00224.10 


SXPhe 


19.353639 


37.710163 


21.72 0.71 


1.06 


s 


59.801 


0.12261 


0.20 


03653_3 


RRc 


19.347147 


37.992413 


17.21 0.57 


0.58 


k 


59.937 


0.32654 


0.39 


00345 _1 


RRab 


19.325082 


37.964170 


18.28 






60.151 


0.57866 


0.72 










Irregular variables 












V92 


IRR 


19.350754 


37.766876 


18.10 0.91 




k 






0.10 


V83 


IRR 


19.346220 


37.737232 


19.10 1.02 


1.05 


k 






0.07 


V93 


IRR 


19.351452 


37.785687 


18.12 0.98 


1.03 


s 






0.04 










Cataclysmic variables 












V15(=B7) 


CV 


19.352057 


37.799019 


18.26 0.20 




k 






0.06 


B8 


CV 


19.343262 


37.747833 


20.64 -0.23 


0.78 


k 






2.27 


06289_9 


CV (?) 


19.348976 


37.770355 


22.80 0.25 


0.88 


s 






3.10 



Table 4. Coordinates and light curve parameters of the contact binaries (W UMa systems, EW) belonging to NGC 6791. V is the 
brightness at maximum and Tq is the time of the primary minimum. 



Star 


Q'2000 


£2000 


V 


<B-V> 


< V - I > 


Ref. 


T 


Period 


Ampl. 


Notes 








[mag] 


[mag] 


[mag] 




[HJD-2452400] 


[d] 


[mag] 




01441_8 


19.339422 


37.778118 


19.98 


1.38 




k 


63.073 


0.24544 


0.07 


likely memb. 


V118 


19.347500 


37.651222 


17.68 


0.75 


1.01 


s 


59.912 


0.30623 


0.70 


likely memb. 


V5 


19.346258 


37.813354 


17.19 


0.90 


0.95 


k 


60.221 


0.31274 


0.05 


member (98%) 


V3 


19.354380 


37.769349 


18.51 


1.05 


1.06 


k 


59.798 


0.31764 


0.09 


member (78%) 


V4 


19.348396 


37.806652 


17.72 


1.01 




k 


59.591 


0.32568 


0.10 


member (98%) 


V8 


19.341938 


37.865810 


17.81 


0.79 


0.88 


k 


59.896 


0.33406 


0.10 


likely memb. 


VI 17 


19.343433 


37.665848 


17.66 


0.87 


0.90 


k 


59.987 


0.36644 


0.38 


likely memb. 


V7 


19.340271 


37.821892 


17.63 


0.93 


0.86 


k 


59.820 


0.39174 


0.31 


likely memb. 



We note that our period for VI 19 is quite different from that 
given by H05 (0.1133 days); the new period makes this star an 
intermediate case between semi-detached and contact systems. 
Error bars on the periods in Tab. [5] are ~ 10 4 d for P <1.0 d, 
~ 10 -3 d for 1.0 < P <2.0 d and a bit larger for P >3.0 d. 

Some variables show very sharp eclipses and out-of-eclipse 
variability due to different levels of stellar activity (05736_9, 
00645_10, V109, 01393_1, VI 1 and V107; for the period of the 
latter star we prefer the longer of the two values given by M05). 

In many cases we observed one eclipse only and we cannot 
give any value for the period, unless it has been given in the 
previous studies, as for the cluster member V80 (86% on the 
basis of the Cudworth membership probability). We also note 
that the amplitude we observed in V80 is much larger than that 
reported by B03. 

To establish the membership of these eclipsing systems is not 
an easy task, since binary effects should be taken into account 
when considering colors and magnitudes. However, on the ba- 
sis of the distance from the cluster center and their position in 
the CMDs, we can argue that V60, 0246 1_8 (both single-event 



eclipsing binaries), 05736_9, V29 and 00645_10 are very proba- 
ble members. This hint is corroborated by the membership prob- 
abilities for V60 and 02461_8, which are 91% and 88% respec- 
tively. 

The special cases of V9 and B4 deserve attention. V9 is the 
binary closest to the center and its membership probability is 
82%. However, it looks a very evolved object in the CMD; its 
period (3.2 d) and activity (see above) are also more typical for a 
Main Sequence star. Therefore, its membership is very doubtful. 

The Cudworth membership probability for B4 is only 40%, 
but in the CMDs B4 belongs to a little "clump" of very blue 
stars. This location is in agreement with the results of Liebert et 
al. ([1994 ) and therefore B4 is likely a blue extanded horizontal- 
branch star belonging to NGC 6791. The star is classified by 
M02 and M03 (who consider it a non-member) as an eclipsing 
binary, but we note that the light curve could also result from a 
rotational modulation. 

Other possible members are: V107, 00331J, V109 and VI 1, 
considering that they are within 6'4 radius from the cluster cen- 
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Table 5. Eclipsing variables with well defined light curves. EA stands for a B Per system, EB for a B Lyr one. V is the brightness at 
maximum and Tq is the time of the primary minimum. Also B4, V60 and 0246 1_8, whose parameters are reported in the Appendix, 
are possible cluster members. The last column shows the Cudworth membership probability (when available). 



Star 


T 

type 


Q '20()0 


x 

02000 


\i 


< B — V > 


< V — I > 


Ref. 


7",, 

1 


Period 


Ampl 


Notes 










[mag] 


[mag] 


[mag] 




THTD-24524001 


L U J 


[mag] 




V119 


EB 


19351961 


37.916328 


18.15 


1.13 


1.33 


k 


59.879 


0.30197 


0.15 


member ? 


V29 


EB 


19354796 


") "7 TCI 10^ 

37.751386 


20.00 


1.23 


1.61 


k 


69.012 


0.43662 


0.22 


likely memb. 


01558_5 


EB 


19.372697 


37.953469 


19.15 








69.028 


0.52910 


0.28 


likely non-memb. 


01393.1 


EA 


19.329588 


37.970392 


21.23 








59.326 


0.58998 


0.56 


likely non-memb. 


0033 1_3 


EB 


19.352367 


37.864391 


19.73 


1.20 


1.36 


k 


68.815 


0.7347 


0.13 


member ? 


Vll 


EA 


19.342575 


37.804802 


19.38 


0.96 


1.22 


k 


67.875 


0.8833 


0.48 


member ? 


05736_9 


EA 


19.348484 


37.721855 


20.20 


1.21 




k 


68.333 


1.210 


0.29 


likely memb. 


V12 


EB 


19.345259 


37.849083 


17.52 


0.96 




k 


64.103 


1.524 


0.06 


member (96%) 


00645_10 


EA 


19.354692 


37.710104 


20.60 


1.33 


1.46 


k 


60.893 


1.451 


0.20 


likely memb. 


V9 


EB 


19.346634 


37.777035 


17.15 


1.23 


1.38 


k 


63.873: 


3.2 


>0.2 


member (82%) 


V107 


EA 


19.355068 


37.761553 


17.97 


0.93 


1.00 


k 


64.433 


3.27 


0.24 


member ? 


V109 


EA 


19.342716 


37.793961 


20.73 


1.46 


1.60 


k 


69.021 


3.70 


0.86 


member ? 



ter. The location in the CMDs of the eclipsing binaries belonging 
to NGC 6791 is shown in Fig. 10 (triangles). 



4.6. Rotational variables 

We found 89 variables whose light curves are characterized by 
small amplitude (usually less than 0.10 mag) and continuous 
variability. It is difficult to ascribe such variability to contact bi- 
naries undergoing grazing eclipses, since they should be less nu- 
merous than those having partial eclipses, since grazing eclipses 
occur only for a particular orientation of the orbital plane. Our 
hypothesis is that in most cases this variability results from spots 
carried by the stellar rotation; under this hypothesis, a large va- 
riety of light curves can be produced. Of course, we cannot rule 
out that a small fraction of these light curves might be actually 
generated by grazing eclipses. 

The complete list of the rotational variables and their light 
curves is given in the Appendix. Here we discuss some exam- 
ples. If the inclination of the rotational axis causes the progres- 
sive disappearance of the largest spots, the light curve displays 
continuous variation, which could be sine shaped in the simplest 
cases (a fraction of the spots is always visible; it can also produce 
Cepheid-like variability, as in the 001606_1 case, Fig. 12 1, or 
with a standstill (the hot or cold spots totally disappear; 00513-2 
in Fig. 12 1 or, more commonly, it can be distorted by other spots 



besides the largest ones (0047 02 in Fig. 12 1. In cases of very 
active stars, a secondary wave also occurs (01 175_5 in Fig. 12 1. 
Since the second wave often covers less than half of the period, 
these rotational variables can be distinguished from eclipsing bi- 
naries; we also note that the amplitude ratio between the first and 
second waves can be very different. 

Also in the three cases in which the full amplitude is larger 
than 0.10 mag (V2, 02006_1 and 07483_9) rotational effects ex- 
plain the observed features better than eclipses. For example, 
the light curves of V2 (P=0.273 d) and 01298_5 (P=0.586 d; see 
Fig. [12 1 show typical eclipsing binary behaviour, but the ampli- 
tudes, the periods, and, mainly, the asymmetries are more typical 
of a rotational effect. The case of 02270_1 1 is different (Fig. 12 1. 



Its light curve is very similar to that of a contact binary, but it 
does not repeat exactly, and unusual scatter is observed through 
the cycle. We also note that this non-repetitive behaviour of the 
light curves, due to the spot activity, is the reason why several 



variables stars show residual standard deviations higher than ex- 
pected. 

Our periods for V34, V37 and V38 are approximately half 
of those given by M02, since these authors classified these vari- 
ables as ellipsoidal ones; the large amplitudes (0.18, 0.06 and 
0. 13 mag) are more in favour of a variability resulting from large 
spots, rather than the purely geometrical effect of tidally dis- 
torted stars. We also note that V37 did not show any flare activity 
similar to that reported by M02 during our survey. We have also 
revised the classification of V16, considered an eclipsing binary 
byM02andM03. 

We count 33 rotational variables in the lO'-circle 
(i.e., 0.105 star/arcmin 2 ) centered on the cluster, while we 
have 56 variables in the remaining 924— arcmin 2 area (i.e., 
0.061 star/arcmin 2 ). We have color indices (B - V and/or V — T) 
for 48 stars; 33 of them have a radial distance less than 10' from 
the cluster center. We can confirm the membership for 6 stars 
having proper motion membership: V16, V38, V42, V48, V53 
and 03079 _9. We have no photometric indices for V41; however, 
it is at only 2' from the cluster center and its Cudworth probabil- 
ity membership is 77%. Therefore, we consider V41 a member. 
For V14 we have the opposite situation because this star is at 
1' from the cluster center and its positions in the CMDs agree 
very well with a membership, but the proper motion measure- 
ments rule out that it can be a cluster member (0%) (see Figure 
13 V14 is displayed as a starred dot). As mentioned by M03, 



the positions of V17 in the CMDs are unusual. Other variables 
located below the subgiant branch like V17 were found in the 
open cluster M67 (Mathieu et al.,|2003) and in the globular clus- 
ter 47 Tuc (Albrow et al.. l20031 l. Probably these objects (named 
"red stragglers" or "sub-subgiant branch stars") are the result of 
some kind of mass exchange between the members of a binary 
system. 

Putting the rotational variables without proper motion mea- 
surements on the CMDs we could infer that 8 stars are lo- 
cated on or close to the MS (represented with filled circles in 
Fig. [T3j ); thus we suggest that these 8 stars belong to the clus- 
ter as well. Among the variables at greater distances, for three 
stars (01149^, 01122_4 and 00513^, all located between 11' 
and 13') the membership is doubtful, since their position in the 
CMDs is unclear. The other stars show apparent magnitudes 
and/or color indices too discrepant to be considered active MS 
stars belonging to NGC 6791. 
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Fig. 7. Top: Light curves of B8 (CV star) and 06289_9 (candidate CV). 
Bottom: positions in the two-colour diagram for two irregular stars 
(V83 and V93), for the cataclysmic variable B8 when in a low state, 
and for the new variable candidate 06289-9. 



When considering the variables without color indices, only 
two (V41 and 01874J2) are at less than 10' from the cluster cen- 
ter. We know that V41 is a probable cluster member (member- 
ship probability 77%), but, at the moment, we have no valid rea- 
son to consider the other star as a member. 

Table [6] lists the rotational variables we suggest as cluster 
members. The error bars on the period are ~ 10 4 d for P < 1 .0 d, 
~ 10 3 d for 1.0< P <2.0 d, ~ 10~ 2 d for 2.0< P <5.0 d; periods 
longer than 5.0 d are tentative. Figure [13] shows the CMDs with 
the rotational variables belonging to the cluster (Tab. [6]i clearly 
indicated. We rejected as cluster members 16 stars out of 32 
located within 10' from the cluster center; i.e., we considered 
them to be stars of the Galactic field. We note that the resulting 
density of the Galactic field (0.051 star/arcmin 2 ) superimposed 
on the cluster is in good agreement with that of the surround- 
ing galactic disk field (0.061 star/arcmin 2 , see above), especially 
considering that Poisson statistics supply uncertainties around 
+0.01 on the density values. 

The stellar rotation and the activity level are both expected 
to be small for single stars as old as NGC 6791. Therefore we 
suggest that the rotational variables belonging to the cluster are 
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Fig. 8. The light curves of the contact binaries that are likely mem- 
bers of NGC 6791. The case of 01434_3 (amplitude much larger than 
0.75 mag) is also shown in the last panel. 



likely short-period binaries, whose rotational velocity and activ- 
ity level have been enhanced by the tidal synchronization. 



4. 7. Long-period variables 

We detected numerous stars having different mean magnitudes 
on the different nights. Their behaviours are more diversified 
than those of the stars we considered as spurious on the basis 
of their close similarities. The resulting power spectra are domi- 
nated by terms at very low frequencies, corresponding to periods 
often much longer than 10 d. These periods cannot be evaluated 
in a precise way, being comparable or, more frequently, longer 
than our time baseline. Therefore, we can only argue that these 
stars are variables, either in a periodic or in an irregular way. 
Since we detected many spotted stars, it is quite obvious to think 
that most of these long-period variables are spotted stars having 
a rotational periods longer than 10 d. The mean amplitude of 
these stars is about 0.02 mag, except for 5 stars whose amplitude 
exceeds 0. 1 mag. 

Among the long-period variables, we used the Cudworth 
probabilities to establish the membership of 18 stars. In order to 
roughly estimate the membership of the remaining long-period 
variables we checked their locations in the CMDs, in the cases 
where at least one color is available. We suggest that 5 stars are 
likely members of NGC 6791: 02138_8, 01610_9, 04392_3, V75 
and 02268_10 (see Figure 14 1. They lie along the MS or the red- 
giant branch and, furthermore, they are all located at distances 
smaller than 8'5, from the cluster center. Looking at the position 
of the variable V76 = V85 (memberships: 97%) in both CMDs, 
we suggest that this star could be similar to the "sub-subgiant 
branch" star VI 7. In Figure [l5]we show its light curve and those 
of the 5 stars that we suspect to belong to the cluster. Table [7] 
lists the long-period variables we suggest as cluster members; 
the entire sample is listed in the Appendix. 
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Table 6. Rotational variables belonging to the cluster. V is the mean brightness value. Tq indicates the time of the maximum 
brightness. The last column shows the membership probability (when available) or our photometric membership. The two uncertain 
cases (V41 and V14) are also listed (see text). 



Star 


Tvoe 


zUUU 




V 


< B - V > 


< V-I> 


Ref. 


T 


Period 


Ampl. 


Notes 


















[HJD-2452400] 








04803.9 


ROl 


19.347698 


37.796043 


21.81 


1.33 


1.84 


s 


61.799 


1.1034 


0.17 


member (B06) 


V82 


ROl 


19.344366 


37.793381 


19.01 


1.00 


1.02 


k 


56.481 


1.1568 


0.04 


likely memb. 


06553.9 


ROl 


19.349134 


37.672577 


19.26 


1.04 


1.13 


s 


61.421 


1.3485 


0.08 


likely memb. 


01724_9 


ROl 


19.344957 


37.785362 


20.73 


1.29 


1.69 


s 


64.410 


1.6130 


0.17 


likely memb. 


V38 


ROl 


19.351021 


37.768288 


18.82 


0.96 




k 


55.630 


1.96 


0.13 


member (92%) 


03079_9 


ROl 


19.346190 


37.754753 


19.23 


1.14 


1.26 


k 


66.630 


2.640 


0.07 


member (93%) 


V14 


ROl 


19.347687 


37.756874 


18.58 


0.93 


1.05 


k 


55.933 


5.45 


0.05 


non— member (0%) 


V48 


ROl 


19.352076 


37.718506 


17.51 


0.88 




k 


65.223 


5.65 


0.09 


member (96%) 


V17 


ROl 


19.344135 


37.817928 


17.92 


1.20 


1.28 


k 


63.211 


6.523 


0.04 


member (88%) 


V51 


ROl 


19.353382 


37.748795 


19.94 


1.22 


1.21 


k 


63.624 


6.72 


0.09 


likely memb. 


V52 


ROl 


19.355795 


37.771935 


17.49 


0.88 


0.88 


k 


64.345 


7.06 


0.03 


likely memb. 


V53 


ROl 


19.350233 


37.743187 


18.72 


0.89 


0.93 


k 


69.294 


7.47 


0.04 


member (86%) 


00436-3 


R02 


19.352205 


37.878635 


18.92 


0.92 


1.08 


k 


60.018 


0.26601 


0.04 


likely memb. 


07483_9 


R02 


19.349997 


37.746311 


21.28 


1.32 


1.70 


s 


60.465 


0.4375 


0.17 


likely memb. 


V41 


R02 


19.347492 


37.806892 


19.09 








60.000 


0.4798 


0.07 


member (77%) 


V42 


R02 


19.350058 


37.714867 


19.51 


1.05 


1.16 


k 


60.323 


0.5068 


0.10 


member (92%) 


V16 


R02 


19.352108 


37.802662 


17.79 


0.93 


1.01 


k 


67.713 


2.182 


0.03 


member (96%) 



Table 7. Long period variable stars that are likely members of NGC 6791 (ordered by increasing right ascension). The column 9 
shows the membership probability (when available) or our photometric membership. 



Star 




#2000 


V 


< B- V > 


<V-I> 


Ref. 


Ampl. 


Notes 


Distance 








[mag] 


[mag] 


[mag] 




[mag] 




arcmin 


02138_8 


19.341864 


37.749653 


19.68 


1.12 


1.18 


k 


0.06 


likely memb. 


4.6 


02444_8 


19.342766 


37.810604 


18.34 


0.93 


1.00 


k 


0.02 


member (78%) 


4.4 


00510_9 


19.343704 


37.796719 


18.67 


0.94 


0.98 


k 


0.03 


member (83%) 


3.4 


01610_9 


19.344818 


37.740486 


19.09 


1.01 


1.05 


k 


0.02 


likely memb. 


3.0 


V94 


19.345139 


37.743549 


17.54 


0.88 


0.92 


k 


0.03 


member (90%) 


2.7 


V95 


19.345295 


37.792412 


19.16 


1.03 


1.10 


k 


0.07 


member (93%) 


2.3 


V56(=V96) 


19.345908 


37.763525 


17.01 


0.95 


0.97 


k 


0.04 


member (98%) 


1.6 


04392-3 


19.345982 


37.870571 


17.90 


0.88 


0.92 


k 


0.02 


likely memb. 


6.1 


V75 


19.346651 


37.766308 


17.38 


0.94 


0.98 


s 


0.01 


likely memb. 


1.0 


04133_9 


19.347109 


37.777020 


18.47 


0.94 


0.98 


k 


0.04 


member (98%) 


0.7 


V76(=V85) 


19.347192 


37.764169 


18.19 


1.03 


1.15 


k 


0.11 


member (97%) 


0.8 


V86 


19.347258 


37.808815 


19.44 


1.06 


1.14 


k 


0.03 


member (83%) 


2.3 


V87 


19.347996 


37.749668 


18.12 


0.91 


0.93 


k 


0.03 


member (98%) 


1.3 


05740_9 


19.348534 


37.808022 


17.96 


0.93 




k 


0.03 


member (95%) 


2.2 


06725 _9 


19.349329 


37.724495 


17.77 


0.91 




k 


0.02 


member (91%) 


3.0 


06796_9 


19.349415 


37.769268 


18.46 


0.92 


1.04 


k 


0.03 


member (92%) 


1.0 


V90 


19.349686 


37.746449 


18.11 


0.86 


0.94 


k 


0.01 


member (94%) 


1.9 


07680-9 


19.350193 


37.718224 


18.32 


0.89 


0.97 


k 


0.03 


member (98%) 


3.5 


V31 


19.350683 


37.785912 


17.12 


1.00 


1.02 


k 


0.01 


member (97%) 


2.1 


09376-9 


19.351952 


37.831886 


18.21 


1.00 


1.02 


k 


0.01 


member (92%) 


4.6 


09611_9 


19.352139 


37.773365 


17.97 


0.85 


0.91 


k 


0.01 


member (76%) 


2.9 


V58 


19.354042 


37.801240 


17.52 


0.89 


0.94 


k 


0.05 


member (87%) 


4.6 


02268_10 


19.359475 


37.731544 


18.59 


0.95 


0.99 


k 


0.02 


likely memb. 


8.5 



5. Conclusions 



Our wide-field survey of NGC 6791 for the planetary-transit 
search allowed us to discover 260 new variable stars. When con- 
sidering the membership probabilities given by Cudworth and 
B06, 13 of them belong to the cluster and one star (09831-9) is 
not member. On the basis of the distances from the cluster center 
and the positions in the CMDs, we suggest that another 1 1 stars 
are likely members, for 22 stars the membership is doubtful, and 



137 stars are likely non-members. No photometric or kinematic 
data are available for 76 stars. 

The variable star content of the cluster is very similar to that 
of the surrounding Galactic environment: in both samples we 
find rotational variables, contact and eclipsing systems. Contact 
binaries and rotational variables belonging to the cluster have the 
same characteristics as those located in the surrounding Galactic 
field. No evidence of pulsating variables has been found in 
NGC 6791, but this is not surprising, since it is a very evolved 
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5 10 15 

radial distance (arcmin) 

Fig. 9. The m — M values calculated for contact binaries by means of 
the Rucinski ( 2003) P-L-C relation plotted against the distance from the 
cluster center. Triangles indicate the stars whose membership has been 
proposed by M03, the filled circles the stars whose membership has 
been proposed by us, and the open circles the stars that we suggest do 
not belong to NGC 6791. The starred point indicates the star 0989 1_9 
which does not belong to the cluster on the basis of the proper motion 
(B06). The (m - M) v value for the cluster (solid line) with an error bar 
of ±0.20 mag (dashed lines) is also shown. 



cluster and stars located in the instability strip or hotter pulsators 
have already left the MS. 

The discovery of the new cataclysmic variable 06289_9 in 
addition to B8 and V15 adds another peculiarity to NGC 6791, 
making it unusual among the open clusters. 
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Fig. 10. The CMDs for the binary systems belonging to NGC 6791. 
Filled circles are contact binaries (W UMa stars), triangles are detached 
or semi-detached systems (J3 Per or f} Lyr systems). 
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Fig. 11. The light curves of short-period detached or semi-detached 
eclipsing binaries. 
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Fig. 13. (V-I)- V and (S- V)-V diagrams for NGC 6791. The rota- 
tional variables that we suggest may belong to the cluster are indicated 
with filled circles. Triangles: stars belonging to the cluster according to 
the Cudworth's membership; starred point: V14, open square: V17 (see 
text for details about these stars). 
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Fig. 12. The light curves of a small sample of rotational variables, il- 
lustrating the growing importance of the second wave. 
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Fig. 14. (6 - V) - V and (V - I) - V diagrams for NGC 6791. Filled circles: long-period variables that we suggest belong to the cluster. Filled 
triangles: long-period variables that belong to the cluster (membership probability higher than 76%.) 
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Appendix A: List of identified variables 

This Appendix includes the full list of the identified variables, 
separated according to our classification: 



L. 
2. 

3. 
4. 

5. 
6. 



Pulsating, irregular and cataclysmic variables: Table: A.l 



EW-Type stars: Table|A.2| Figure A.l 



EA/EB-type stars: Table|A.3 

Single-wave rotational variables: Tables |A.4| |A.5 



Figures A.2 A.3 



Double-wave rotational variables: Table A. 6 Figure A.4 



Long-period variables: Table A. 7 



Into the tables, for each star we give the name (a five-digit 
number followed by the chip number which the star belongs), 
coordinates, photometric data (always the V mag, B — V color 
when available), informations about the variability (To, period, 
amplitude), distance from the center (in arcmin) and finally the 
numerical value of the Cudworth's membership probability (re- 
ported in the column 'Memb.'). 

In most cases, when membership probabilities were not 
available, in the same column the label "m" means that we re- 
tain the star belonging to the cluster, while 'm?' and "nm" mean 
"uncertain membership" and "likely non-member" respectively. 
The label "ndl" means that no photometric data are available to 
advance hypothesis about the membership, but the star is located 
nearer than 10' from the center of the cluster. Finally "nd2" 
means that no photometric data are available and the star is lo- 
cated further than 10' from the center; in this case we strongly 
suggest that the star does not belongs to the cluster. 
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Table A.l. Pulsating, irregular and cataclysmic variables. V is the minimum brightness for CVs and irregular, the mean brightness 
for pulsating variables. r () is the time of maximum brightness for pulsating stars. 



Star 


Type 




£2000 


V 


< B - V > 


< V - 1 > 


Ref. 


T 


Period 


Ampl. 


Memb. 


Distance 










[mag] 


[mag] 


[mag] 




[HJD-2452400] 


[d] 


[mag] 




[arcmin] 


V123 


HADS 


19.362064 


37.666034 


17.08 


0.45 




k 


59.559 


0.06026 


0.14 


nm 


11.8 


01497.12 


HADS 


19.379083 


37.812419 


16.06 








59.528 


0.07227 


0.40 


nd2 


22.2 


00311_7 


SXPhe 


19.324628 


37.716768 


23.17 








59.605 


0.10443 


0.10 


nd2 


17.0 


00224.10 


SXPhe 


19.353639 


37.710163 


21.72 


0.71 


1.06 


s 


59.801 


0.12261 


0.20 


nm 


5.4 


03653_3 


RRc 


19.347147 


37.992413 


17.21 


0.57 


0.58 


k 


59.937 


0.32654 


0.39 


nm 


13.3 


00345_1 


RRab 


19.325082 


37.964170 


18.28 








60.151 


0.57866 


0.72 


nd2 


20.0 


V92 


IRR 


19.350754 


37.766876 


18.10 


0.91 




k 






0.10 


m 


1.9 


V83 


IRR 


19.346220 


37.737232 


19.10 


1.02 


1.05 


k 






0.07 


m 


2.4 


V93 


IRR 


19.351452 


37.785687 


18.12 


0.98 


1.03 


s 






0.04 


m 


2.6 


V15(=B7) 


CV 


19.352057 


37.799019 


18.26 


0.20 




k 






0.06 


98 


3.3 


B8 


cv 


19.343262 


37.747833 


20.64 


-0.23 


0.78 


k 






2.27 


m 


3.7 


06289.9 


CV 


19.348976 


37.770355 


22.80 


0.25 


0.88 


s 






3.10 


m (B06) 


0.7 



Table A.2. Contact binaries; V is the brightness at maximum and Tq is the time of the primary minimum. 



Star 


^2000 


^2000 


V 


< B- V > 


< V- I > 


Ref. 


T 


Period 


Ampl. 


Memb. 


Distance 








[mag] 


1 mag ] 


1 mag ] 




[HJD-2452400] 


[d] 


[mag] 




[arcmin] 


V122 


19.360729 


37.641436 


20.92 








59.614 


0.22883 


0.58 


nd2 


11.9 


V115 


19.330646 


37.975902 


20.70 








59.473 


0.23636 


0.24 


nd2 


17.4 


01407-8 


19.339256 


37.701576 


21.52 


0.64 


1.25 


k 


59.447 


0.24155 


0.82 


nm 


7.5 


01150.4 


19.357227 


37.962177 


18.58 


1.07 


0.93 


k 


59.455 


0.24510 


0.27 


nm 


13.2 


01441.8 


19.339422 


37.778118 


19.98 


1.38 




k 


63.073 


0.24544 


0.07 


m 


6.2 


00144-2 


19.334044 


38.030701 


19.19 








59.457 


0.24780 


0.77 


nd2 


18.5 


01670.10 


19.357625 


37.681442 


16.64 


1.14 


1.24 


k 


59.729 


0.25807 


0.49 


nm 


8.7 


V121 


19.358063 


37.932346 


17.24 


0.81 


0.84 


k 


59.619 


0.26742 


0.71 


nm 


12.0 


0229 111 


19.371624 


37.795464 


19.15 








59.582 


0.26774 


0.29 


nd2 


16.8 


V23 


19.338614 


37.787781 


16.92 


1.04 


1.23 


k 


59.915 


0.27180 


0.07 


nm 


6.8 


V25 


19.328426 


37.713237 


18.50 








59.772 


0.27730 


0.56 


nd2 


14.4 


00665_12 


19.375697 


37.713244 


18.89 








59.472 


0.28369 


0.41 


nd2 


20.0 


09831_9 


19.352356 


37.780907 


20.55 


1.09 


1.19 


s 


59.662 


0.29488 


0.3 


nm (B06) 


3.1 


01434.3 


19.350643 


37.985512 


22.10 








59.831 


0.30581 


0.97 


nd2 


13.0 


V118 


19.347500 


37.651222 


17.68 


0.75 


1.01 


s 


59.912 
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Table A.3. Eclipsing variables. V is the brightness at maximum and 7q is the time of the primary minimum. 
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Table A.4. Rotational variables with a single-wave light curve. V is the mean brightness value. Tq indicates the time of the maximum 
brightness. 
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Table A.5. Continue from Table 



A.4 



other rotational variables with a single-wave light curve. V is the mean brightness value. Tq 



indicates (one of) the times of maximum brightness. 



Star 




•feooo 
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< B -V > 


< V - I > 
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Period 
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Memb. 


Distance 
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[mag] 




rii 1 1"\ A r" r\ A i~\f\~\ 
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Table A.6. Rotational variables with a double-wave light curve. V is the mean brightness value. 7q indicates (one of) the times of 
maximum brightness. 
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[arcmin] 


00436.3 


19.352205 


37.878635 


18.92 


0.92 


1.08 


k 


60.018 


0.26601 


0.04 


m 


7.1 


V2 


19.354875 


37.766689 


19.74 


0.93 


1.21 


k 


59.712 


0.27344 


0.17 


nm 


4.9 


02006.1 


19.332258 


38.043080 


21.54 








60.448 


0.37500 


0.19 


nd2 


19.8 


01175.5 


19.370209 


37.998342 


19.73 








59.918 


0.41481 


0.09 


nd2 


20.8 


00536.11 


19.365067 


37.716529 


19.13 








60.460 


0.43740 


0.03 


nd2 


12.6 


07483.9 


19.349997 


37.746311 


21.28 


1.32 


1.70 


s 


60.465 


0.4375 


0.17 


in 


2.1 


V41 


19.347492 


37.806892 


19.09 








60.000 


0.4798 


0.07 


77 


2.2 


V42 


19.350058 


37.714867 


19.51 


1.05 


1.16 


k 


60.323 


0.5068 


0.10 


92 


3.7 


01362.7 


19.329000 


37.662097 


20.53 








60.502 


0.5560 


0.08 


nd2 


15.1 


02270.11 


19.371548 


37.794338 


17.95 








60.525 


0.5679 


0.05 


nd2 


16.8 


01298.5 


19.371044 


38.002255 


19.76 








60.402 


0.5859 


0.06 


nd2 


21.4 


V16 


19.352108 


37.802662 


17.79 


0.93 


1.01 


k 


67.713 


2.182 


0.03 


96 


3.4 


00568.2 


19.336329 


37.881542 


18.52 








60.007 


2.704 


0.08 


nd2 


10.6 


00019.7 


19.323382 


37.710172 


18.95 








65.473 


5.882 


0.12 


nd2 


17.9 



De Marchi et al.: Variable stars in the open cluster NGC 6791 and its surrounding field 



P = Q,3l d 

17.i& 




024)3.3 P = 0.38 d 00513.S P = Q,45 d 

Si,6 r ; i 



02393.10 P=0-48 d 0] 487.11 F=G.S? d 



a ao.is ! 



A 



0177«_4 P-D.6B d 

I9.T? 




] O Q.5 I 

Q03OI 5 P-0.7D d OODBH_H P-D.ll d 

3S 




D.5 1 
01IO6 IB P-O.BE d ODB13_B P-l.OB d 




23.1 ■: W\ "4 



■ t 151 
7 BB 



r 1 — i ■ ■ — ■ ■ I i — i — i i — |- 1 - : 



A 



2D. 9 



□484)3-9 P-l 10 d OlBOe.ll P-I.IZ d 

so.: 



i cu t - 

v - y 



(J U_!i 

vaa p=i.ia d 

i«.as K 1 " "j : ' 



□ d.& l CO!,! 
V34 P = l 20 d O5302_3 P=1.33 d 




0.3 k J U 

_JE3 so.4s £-■"> .... i ... j 



19 



1&.06 ; 







(f 6 1 
00110.5 P=1.41 d 




0.6 ] 0.6 I 

# if 
066A3 S P=1.36 d 00471_lS P=l 35 d 

IB. 2 b- ' I ' 1 ' ' I - 

19.22 [ /jL -I BO. IE | 

SB.-24 I *i 4 80.14 E- 

T9.2B [ / 30.iai 

■j ■ !: v« -^jg ao.is I 
■ Tjr V 3*0.2! 



CO ] W 0.6 1 

Ol074_3 P = 1.SS d VI 11 P=l.Sfl d 

■t^ 20.55 

20. a 

20-$* i 




* ,...jH 



0.0 

vst p=i.ai d 




0.6 

0l?S4_9 P=l.ei d 



(J Q-Si 1 D.S ] 0.5 1 

£..10 P-l-87 d 011W.11 P- 1 fZ d OQQ1H.5 P- 1 79 a 




19.26 
19.2a 

1 B.3 
19.32 
I9.3-5 

i9.aa 




0.6 ] 0.6 

00738.13 P- 1 ,015 d 00771-11 F-1.fl5d 

-i — | — i — r 

:9 32 
]93+ 
■9 36 



? is- 

/ ■ 



is. s r 



1 1M y\ 



0.5 

00576. IS P-S,31 

"9 By r 1 | 1 
19.64 I i 

■ i i.u --j- 

■J !i f 

BB 



0.& 



an 



- 



j — | — i — i I i i 



D.S 



Fig. A.2. Rotational variables with a single-wave light curve, in some cases very distorted. 
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Fig. A.3. Rotational variables with a single-wave light curve, in some cases very distorted. 
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Fig. A.4. Rotational variables with a single-wave (from V51 to 05877_9) and a double-wave light curve. 
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Table A.7. Long period variable stars (ordered by increasing right ascension). 
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